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Abstract. The use of levitated nanospheres represents a new paradigm for the 
optomechanical cooling of a small mechanical oscillator, with the prospect of realising 
quantum oscillators with unprecedentedly high quality factors. We investigate the 
dynamics of this system, especially in the so-called self-trapping regimes, where one 
or more optical fields simultaneously trap and cool the mechanical oscillator. The 
determining characteristic of this regime is that both the mechanical frequency ojm 
and single-photon optomechanical coupling strength parameters g are a function of 
the optical field intensities, in contrast to usual set-ups where ujm and g are constant 
for the given system. We also measure the characteristic transverse and axial trapping 
frequencies of different sized silica nanospheres in a simple optical standing wave 
potential, for spheres of radii r = 20 — 500 nm, illustrating a protocol for loading single 
nanospheres into a standing wave optical trap that would be formed by an optical 
cavity. We use this data to confirm the dependence of the effective optomechanical 
coupling strength on sphere radius for levitated nanospheres in an optical cavity and 
discuss the prospects for reaching regimes of strong light-matter coupling. Theoretical 
semiclassical and quantum displacement noise spectra show that for larger nanospheres 
with r > 100 nm a range of interesting and novel dynamical regimes can be 
accessed. These include simultaneous hybridization of the two optical modes with the 
mechanical modes and parameter regimes where the system is bistable. We show that 
here, in contrast to typical single-optical mode optomechanical systems, bistabilitics 
are independent of intracavity intensity and can occur for very weak laser driving 
amplitudes. 



PACS numbers: 03.75.Lm, 67.85.Hj, 03.75.-b, 05.60. Gg 
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Introduction 

Extraordinary progress has been made in the last half-dozen years [H [2] towards the goal 
of cooling a small mechanical resonator down to its quantum ground state and hence 
to realise quantum behavior in a macroscopic system. Implementations include cavity 
cooling of micromirrors on cantilevers [3j HI EJ E] ; dielectric membranes in Fabry- Perot 
cavities [7]; radial and whispering gallery modes of optical microcavities [8] and nano- 
electromechanical systems j9]. Indeed the realizations span 12 orders of magnitude in 
size j2], up to and including the LIGO gravity wave experiments. In 2011 two separate 
experiments [TOj, [H] achieved sideband cooling of micromechanical and nanomechanical 
oscillators to the quantum ground state. In Ref. [12], spectral signatures (in the form of 
asymmetric displacement noise spectra) of quantum ground state cooling were further 
investigated. Corresponding advances in the theory of optomechanical cooling have also 
been made p31 [151 UM E] • 

Over the last year or so, a promising new paradigm has been attracting much 
interest: several groups [HI [191 1201 I2H (221 [23] have now investigated schemes 
for optomechanical cooling of levitated dielectric particles, including nanospheres, 
microspheres and even viruses. The important advantage is the elimination of the 
mechanical support, a dominant source of environmental noise which can heat and 
decohere the system. 

In general, these proposals involve two fields, one for trapping and one for cooling. 
This may involve an optical cavity mode plus a separate trap, or two optical cavity 
modes: the so-called "self-trapping" scenario. 

Mechanical oscillators in the self-trapping regime differ from other optomechanically- 
cooled devices in a second fundamental respect (in addition to the absence of mechanical 
support): the mechanical frequency, associated with centre-of-mass oscillations is 
not an intrinsic feature of the resonator but is determined by the optical field. In partic- 
ular, it is a function of one or both of the detuning frequencies, A x and A 2 , of the optical 
modes. In previous work [23], we analysed cooling in the self-trapped regime and found 
that the optimal condition for cooling occurs where both fields competitively cool and 
trap the nanosphere. This happens when ojm is resonantly red detuned from both the 
detuning frequencies i.e. um(^i, A 2 ) ~ — Ai j2 so the relevant resonant frequencies are 
mutually interdependent. Most significantly, the effective light-matter coupling strength 
g also depends on the detunings. 

The effective coupling strength, g = g^fri (the optomechanical coupling rescaled 
by the square root of photon number) determines whether one can attain strong 
coupling regimes in levitated systems such as recently observed in a non-levitated set- 
up [21]. It determines too whether one may access other interesting dynamics, both 
in the semiclassical and quantum regimes. We consider in particular the possibility of 
simultaneous hybridization of the two optical modes with the mechanical mode; here, 
we consider also the implications of a static bistability, which, unusually, occurs also in 
the limit of quite weak driving in the levitated self-trapped system. 
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In the present work, we investigate theoretically and experimentally the strength of 
the optomechanical coupling. In particular, we present experimental measurements of 
the mechanical frequency of a nanosphere trapped in an optical standing wave in order 
to investigate the optical coupling as a function of the size of the nanosphere. 

In section 1 we review the theory of the cavity cooling and dynamics of a self-trapped 
system, and in section 2 we employ the experimentally measured size dependence of the 
coupling to determine the range of optomechanical coupling strengths accessible in a 
cavity. The data suggests that the most effective means to attain stronger coupling will 
be to employ larger nanospheres of typical radii r = 200 — 300 nm. Our work suggests 
that increasing photon number by stronger driving (and by implication increasing 
rescaled coupling strengths) will not prove an effective alternative, since in the present 
system we show g oc n 1//4 rather than y/n, so the rescaled coupling increases very slowly 
with laser input power. 

In section 3, we in investigate the cooling and dynamics. In sec. 3.1, we review the 
corresponding cooling rate expressions obtained from quantum perturbation theory (or 
linear response theory). In sec. 3.2 we report a study of the corresponding semiclassical 
Langevin equations and compare them with fully quantum noise spectra; we compare 
also quantum, semiclassical and perturbation theory results for levitated nanospheres. 
In section 4 we investigate novel regimes of triple mode hybridization, coincident with 
static bistabilities, which the present study shows are experimentally accessible given the 
large optomechanical coupling strengths associated with r = 100 — 300 nm nanospheres. 
In section 5 we describe the experimental study which provides data from which the 
size- dependence of the coupling may be inferred. In Section 6, we conclude. 

1. Theory: Quantum Hamiltonian for a nanosphere in a cavity 

We approximate the equivalent cavity model by a one-dimensional system, with centre- 
of-mass motion confined to the axial dimension. In this simplified study, we consider 
only the axial dynamics: for the cavity system, we will have a much smaller tranverse 
frequency relative to the axial frequency, i.e. u t <C uj a and there is little mixing between 
these transverse and axial degrees of freedom. 

We consider the dynamics of the following Hamiltonian: 



Two optical field modes of a high finesse cavity ar e coupled to a nanosphere 
with centre-of-mass position x. The parameter A (dependent on the nanosphere 
polarizability), determines the depth of the optical standing- wave potentials. We 
investigate the case where both modes competitively cool and trap the nanosphere, 
in contrast to previous schemes [HI [19] where one optical field is exclusively responsible 
for trapping, while the other is exclusively responsible for cooling. H is given in the 
rotating frame of the laser which drives the modes with amplitudes E\ and Ei = RE\ 




(1) 
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respectively, where R represents the ratio of driving amplitudes for the two modes. We 
restrict ourselves to the regime < R < 1, since we consider the most general case 
where both optical modes contribute to the trapping as well as the cooling. Thus we 
can define mode 1 simply as the mode which is more strongly driven and mode 2 as the 
mode which is (except where R = 1) more weakly driven. The detunings Aj = oj j l — u J c 
for j = 1,2 are between the input lasers and the corresponding cavity mode of interest, 
and 0i j2 represents the phases of the optical potentials. 

The two fields could represent two modes generated by the same laser field, or they 
could be generated by two independent lasers. Nonetheless, since the particle motion 
is confined to within one wavelegth, one can make the approximation k\ ~ k 2 = k. 
Previous studies generally consider <pi = 0,02 — 7r /4 to be convenient, since then the 
anti-node of one field coincides with a purely linear potential of the other optical field, 
but we may also consider general values of 0i — 02- One can write corresponding 
equations of motion: 



x — y, a j a j sm 2{kx — <pj) — 1 mx 

j 

dj = i^jdj — iEj + iAdj cos 2 (kx — <pj) — —dj, (2) 

where j = 1, 2 for the two optical- mode realisation. Additional damping terms have 
also been added: |oj accounts for photon losses due to mirror imperfections and the 
V mx term for mechanical damping. The above should also include quantum noise terms 
arising from (say) shot noise or gas collisions: for brevity, the quantum noise terms are 
left out until sec. 3. 

We consider here the linearised dynamics; we replace operators by their expectation 
values and perform the shifts about equilibrium values such as dj(t) — > aj + dj(t), and 
x — > xo + x(t). The values for the equilibrium photon fields (e.g. for the two-mode case) 
are «i = —iE 1 [f — iAf] and a 2 = —iRE 1 [f — iAf] . The equilibrium position 
is then given by the relation — = | «2 1 2 / 1 c^i 1 2 ? by numerical solution of the 

equation, 



_ sin2(fcs-fr) = 2 \J-tA^l 

sin2(kx-(j) 2 ) |f -iAf| 2 ' 1 ' 

where AJ = A^ + Acos 2 (kxQ — <pj). As usual we consider the dynamics of the 
fluctuations via the linearised equations. To first order, the linearised equations of 
motion, in the shifted frame, are: 



hkA 



x = -u M x- 



m 



^^(a*dj + ctjd^) sin2(A;:ro — 4>j) — Tm% 



K 



dj = iAjdj — ikAotjX sin2(kx — <pj) — —dj. 



(4) 
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,;.r\, = y, \ a j\ cos2(kxo — 4>j). (5) 



m 



We can restrict ourselves to real equilibrium fields. We take ctj = aje J then 
transform dj — > Oj-e 1 ^. Thus (a^a, + ayS^) = OLjiflj + aj). We also rescale the mechan- 



ical oscillator coordinates x — > \/2X zp fX and p — >■ y/hmuMp, where X 2p _f = y 2m ^ M is 
the zero-point fluctuation length scale. Hence, ^ + \muj 2 M x 2 — > tj ^ 2 M -[x 2 + p 2 )- 

Below we drop the tilde so the equilibrium field values ctj = aj are real. Using field 
operators x = (b + &)/\/2, the linearised dynamics for a two-optical mode system would 
correspond to an effective Hamiltonian: 



^ = - Afdldi - A x 2 d\d 2 + u M (A x , A x pb 
h 

+ <7i(Ai, Af)ai(ai + aI)(S + 6+) + <&(AJ, Af)a 2 (a 2 + a|)(S + 6 f ).(6) 




Figure 1. (Colour online) Size dependent effects in the magnitude of the 
optomechanical coupling parameter g±. It is assumed that cavity parameters would 
correspond to g\ = 10 6 Hz at r = 150 nm, for photon numbers n\ — 10 9 . For 
comparison, the value of A is also shown, as are the experimental and simulated 
frequencies cu a = % = 27t/m- The u>m are scaled by a factor of 10 for clarity and in 
fact, values of ujm ~ 1 MHz are quite realistic in optical cavities. 
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2. Towards strong light-matter coupling with two optical cavity modes 

The Hamiltonian in Eq. [6] appears analogous in form to standard, well-studied 
optomechanical Hamiltonians, albeit with two optical modes rather than one. However, 
it differs in one important respect: in this case, both the mechanical frequency 
u)m(Ai, Af) and the optomechanical coupling strengths gxp = (71,2 (Af, Af ) are not fixed 
and depend on the detunings. The fact that the frequencies of the three modes (two 
optical, one mechanical) are interdependent makes the dynamics different from other 
optomechanical set-ups, where the equilibrium mechanical frequency (i.e. excluding 
shifts arising from the fluctuations) is intrinsic to the mechanical oscillator. 

There is considerable interest in achieving strong-coupling, which leads to regimes 
of light-matter hybridization. The corresponding mode splitting has been observed 
experimentally [24J . In typical set-ups, these regimes are reached if the rescaled effective 
optomechanical coupling exceeds the damping rates i.e. g = g^/n > K,T M , where 
n ~ \a\ 2 is the cavity photon number. Even if the unnormalized coupling is weak, 
strong-coupling regimes may be achieved by increasing the driving power and thus 
increasing intracavity photon numbers. 

In the present levitated system, a particularly interesting regime would involve 
triple-mode hybridization enabling, for example, the coupling of the two modes 
of light via the mechanical mode. However, here, mode hybridization (for which 
9i,2 = (?i,2«i,2 ^ k, T M ) depends non-trivially on the detunings. 

We argue that large coupling cannot be easily achieved by increasing the driving 
power, since gi^ oc n 1 /^ and thus increases slowly with the driving strength. We can 
show, by a simple argument, that increasing the nanosphere size provides the most 
effective means to attain strong coupling. 

For the self-trapped system, optomechanical coupling strengths are gj = 
s/2kAX zp f sin2(fcxo — 4>j) and depend on the detunings via Xq. Note also that 
X zp f = y^h/ (27tium) here too depends on the detunings via %. 

For triple mode hybridization, um ~ Ax ~ A 2 . For convenience, we also take 
0i = 0,02 = 7r /4. Then, since tan2A;x = |a 2 | 2 /l a i| 2 > we can re-write Eq. [5j 



2 *2hAk I .0 

w m = TTi — • (7) 

mcos2fca;o 

We consider near symmetric driving of the two optical modes for which R ~ 1 and 



thus /, ./ u ?a 0/2 = tt/8 so u m ~ (^T") ' n \a - Hence, 



Since the optomechanical coupling increases only very slowly with cavity photon 
number the most effective means to reach strong coupling regimes is to increase the 
nanosphere size to the maximum practical value (r ~ 200 — 250 nm). 
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For the ideal case where the nanosphere radius r is small [TH] ie. A 3> r, then 
A{r <C A) = A (r) where the small nanosphere coupling takes the form: 

where V s = 4/3irr 3 is the sphere volume (and hence m = V s p where the density 
p = 2000 Kgm" 3 for silica). In turn, V c = n(w/2) 2 L is the cavity volume, where 
w ~ 40 fim is the cavity waist and L ~ 0.5 — 1 cm is the cavity length. 

For larger nanospheres, the measured size- dependent corrections must be applied. 
In the experiments described below, we find that the mechanical oscillation frequency 
is modulated by a finite size correction %(r) = u)m{t — 0)/(r) (see Fig. [TTJ and 
description of the measurement of f(r) in section 5 below). Thus, since: 

u4 oc (10) 
m 

then A(r) = A (r)f 2 (r) and the coupling is in turn modulated by the finite size 
correction. The experimental results suggest that for r < 200 nm, then f(r) ~ 1 and 
A(r) ~ A (r). 

For example, for r = 150 nm, L = 1 cm and w = 40 /im, then Aq — 8 x 10 5 Hz. For 
reasonable values of cavity decay constants k~2-8x 10 5 Hz, then for ni ~ 10 9 , 

~ 3l 2 ~ 5.4 x KT 6 (^hd) ~ io 6 Hz > « (11) 
\ m J 

For r < 200 nm, ^ lj2 oc r 3 ' 2 . Thus a 200 nm sphere provides an optomechanical 
coupling about an order of magnitude larger than a 40-50 nm sphere. To achieve 
a comparable increase in coupling by photon number enhancement would require 
increasing the driving power by a factor of order 10 4 . 

A more careful analysis, including the effects of the finite-size correction function 
f 2 (r) is shown in Fig. [TJ We see that g\ (and for m ~ ri2, also g^) reaches a maximum 
value for r ~ 300 nm before falling to zero. Other maxima for larger r do not provide 
a larger value of g\p- Furthermore, they have the disadvantage that they may enhance 
photon recoil heating effects. For comparison, the value of A is also shown, overlaid on 
the experimental and simulated frequencies. 



3. Dynamics 



3.1. Optomechanical damping 

A previous study [23], using rescaled coordinates, investigated the full parameter space 
of two optical mode cooling. Here we investigate more carefully the effect of non-zero 
mechanical damping. Using linear response theory, we can extract cooling rates from 
the equations Eq. HJ 



(12) 
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where 



1 2 2 
\Oij\ 9jK 



SJu) = s-, (13) 

for j = 1, 2. Net cooling occurs for T opt < 0. Although the above is quite similar in 
form to standard optomechanical expressions, as explained previously, rather different 
behaviour is observed since here ujm and gj are both dependent on the AJ. 

From quantum perturbation theory we can show that R n ^ m , the rate of 
transition from state n to n + 1 is: R n ^ n+ i — (n + 1) (Si(um) + S 2 (um)) while 

R n ->n~i = n (Si(-u M ) + S 2 {-com)) 

For n » 1, then R n ^ n+ i — _R n ^ n _i gives the cooling rate of Eq{12j However, with 
the exact expressions we can show that the equilibrium mean phonon number is 

/v = S^um) + S 2 (um) , 14 , 

PT Si(-u M ) + S 2 (-u M ) - Si(u M ) - S 2 (lom)' 




Figure 2. (Colour online) Maps of cooling rates calculated from Eq.[T2]for parameters 
R = 1.0 and R = 0.5. Blue corresponds to cooling, yellow/white to heating. The 
white lines indicate the locus of the single held resonances (where —Af = ljm or where 
— = wjf). The detunings are given in units of A and are dimensionless. For R = 1 it 
is clear that there is a deep, maximum cooling region at a double resonance where the 
two white lines intersect and both optical fields cool simultaneously. It is also evident 
that there is a strong cooling resonance for +Af 2 ± = w m- For R = 0.5, three cooling 
resonances — A^ ± = u>m, — Af-t = and +Af = ujm merge to give a very broad 
strong-cooling region, quite insensitive to detuning A2 over a range of over 1 MHz. 
Here A = k/2 = 0.3 MHz and the input power into mode 1 corresponds to 2mW. 

In Fig.[2]we show colour maps comparing the cooling and minimum phonon numbers 
for both R = 0.5 and R = 1 respectively. The cooling behaviour was investigated 
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previously in [23]. In this case, for each fixed detuning Ai there are up to three cooling 
resonances (at three different values of A 2 ), where strong damping is observed (and 
similarly for each fixed A2). This is in contrast to single optical mode schemes where 
there is a single cooling resonance for which Aj = —u>m or A 2 = For the R = 0.5 

map the three cooling resonances merge, giving a single extended cooling region of about 
1 MHz width. 

For R = 1 the map has a high degree of symmetry, since the role of the two optical 
modes is interchangeable. The figures show that the largest cooling rates are found in 
the double resonance region, making it the most favourable region to work in. 

The equilibrium phonon number in Eq. HH concerns only the idealised situation 
where there is a very good vacuum, negligible photon recoil heating and thus no 
mechanical damping or heating effects. For small r < 200 nm spheres, we assume recoil 
heating is negligible [19] and the dominant source of mechanical damping is background 
gas collisions, which provide an effective mechanical damping T M = ^n^-r 2 n g v g [T8|[T9] 
where m g /m s is ratio of the gas particle's mass to that of the sphere, n g is the gas number 
density and v g is the mean gas velocity for a room temperature thermal distribution. 

It can be shown that the perturbation theory argument above can be adapted to 
obtain equilibrium phonon numbers for a given cooling rate T opt : 



( n )PT = -p 1 ip I , ( 15 ) 

1 M + I 1 opt\ 

where Tb — 300K. Alternatively, the final equilibrium temperatures T eq = 
M r ^_^^ " ac , where T vac is the equilibrium oscillator temperature which would have 
been obtained in a perfect vacuum. 



3.2. Quantum and semiclassical noise spectra 

Although we investigate only a two optical mode system, generalization to more optical 
modes is straightforward. We consider a set of equations of motion, for j = 1, ...N: 

t = -(za; M (Af,...AJ) + ^)6 + ^^(A^...A|)(a j + a]) + V^&m 

j 

(16) 

where the optomechanical strengths ^(A^, ...A|) = — kAX ZPF sin 2(kx — <pj) 
depend on the detunings (as does the mechanical frequency com)- In the two mode 
case we consider, we take 0i = and 2 = 7r /4. 

The optical modes are subject to photon shot noise, while the mechanical modes 
are subject to Brownian noise from collisions with gas molecules in the cavity. For the 
photon shot noise, we assume independent lasers and uncorrelated zero temperature 
noise for which (a\ n (t')a in (t)) = 0, while (a^(t')d^ n (t)) = 5(t — t')5ij. For the gas 
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collisions, we take (b in (t')b\ n {t)) = (n B + l)6(t - t') and (b\ n {t')b in {t)) = n B S(t - t') 
where the number of surrounding bath phonons n B m 

The above equations can be integrated in frequency space to obtain analytical 
expressions for the displacement noise spectra for the arbitrary mode case. We can 
evaluate the displacement spectrum S xx (u) = (\x(oj)\ 2 )qm = ^ J e~ tWT (x(t + T)x(t))dr . 
We obtain: 



(\x(u)\ 2 ) QM \M(u)\ 2 = T M [\ X M(u)\ 2 n B + \x M {-u)\ 2 {n B + 1)] 

+ fiM^)i 2 E^ 2 i^(- w )i 2 ' ( 17 ) 



i=i,2 



where the represent optical and mechanical susceptibilities: 



Xjo{UJ) 



-i(oj + A X A 



K 



-1 



- uj M ) + 



M 



with = Xm{u) ~ X*m(~u) and r]j(uj) = Xjo{u) - X* jo {-u)\ then also 

M{u) = i + ^)Y J3 g 2 3 \nM\ 2 - 



in 



^3 6 

H 

H 4 
CO 

2 




CO 



le-05 



pressure '= 1 mBar 



Quantum spectrum 

Semiclassical spectrum 



pressure = 10 mBar 



Quantum Limit 



J I ^ 



-1 1 

u (2tt x MHz) 



Figure 3. (Colour online) Comparison between quantum and semiclassical 
displacement spectra for gas pressures of 1 mBar and at near vacuum pressure in a 
strong cooling region. At high vacuum, the ground state is approached and thus, for 
the quantum spectrum, the blue sideband vanishes. At higher pressure there is good 
agreement between the quantum and classical results. Spectra near double resonance 
for input power P 1 = 7mW, A = k = 3 x 10 5 Hz, A 1 = —1.5 MHz, A 2 = -0.68 MHz, 
R = 0.5. Some hybridization between the mechanical mode and optical mode 1 is seen 
in the characteristic double-peak sideband structure. 

We compare the quantum displacement with corresponding semiclassical solutions 
in the steady state. The linearised two mode system Eq. El in matrix form corresponds 
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-2 -1 
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Figure 4. For the broad cooling region formed from three overlapping resonances seen 
in Fig. [U[b), we show equilibrium phonon numbers obtained from Eqs. (|15[) . (|17[) and 
(|20p i.e. perturbation theory, the analytical quantum noise formula and semiclassical 
Langevin equations respectively. Agreement between quantum and semiclassical 
results is excellent, less so for perturbation theory at low pressures. 



to a standard problem [30]. Inclusion of the noise arising from gas collisions or laser 
shot noise yields a set of corresponding Langevin equations: = AX + BE(t), where 

A is termed the drift matrix. Its eigenvalues give the stabilities and eigenfrequencies 
of the system's normal modes, while the noise is determined by B, a constant diagonal 
matrix. The elements of the random noise matrix are assumed to be 5-correlated 
(Ei(t)Ej(t')) = 5(t — t')5ij. Methods for obtaining the solution for the steady state 
correlation functions of this system, under conditions of stability, i.e. where all the 
eigenvalues of A have negative real parts, are well-known [30]. The required noise 
spectra, or autocorrelation functions, in frequency space are: 

1 f°° 

S(w) = — / e-*" (X(i + r)X T (i)>dT, (19) 

27T-DO J 

where, 

S(w) = -L (A + luiy 1 BB T (A T - luiy 1 , (20) 

where the diagonal matrix BB T has elements ((rie + \)^m, ( n B + |)Tm, f , f , f , f ) • 
From the above, the noise spectra of all modes may be calculated. Eq.SSC yields 
semiclassical sideband spectra, symmetrical in u. 

In Fig. [3], we compare semiclassical displacement spectra calculated from Eq. [20] 
with corresponding quantum results obtained from Eq. [TT] At high pressures (and hence 
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Figure 5. Triple mode splitting. For A = n = 3 MHz, even at quite high pressures 
(here ImBar), mode splitting is seen in the noise spectra of the optical modes. In all 
the plots, Ai = —1.15 MHz is held fixed while A2 is swept from to -1.6 MHz (for 
an input power of 2mW into mode 1, while R = 0.5). (a) Shows noise spectra for 
both optical mode 1 and mode 2. Three way hybridization between the mechanical 
and both optical modes appears clearly (highlighted in the bold blue line). For clarity, 
some of the strongest peaks have been truncated in height. In (b) three avoided 
crossings are apparent. The dominant character of each normal mode is indicated by 
the colour (black is mechanical, blue is optical mode 1, red is optical mode 2). When 
— A2 is large, there is no mixing. However as A2 — > 0, there is strong mixing and the 
dominant character of each normal mode changes from light to matter (or vice versa) 
as an avoided crossing is encountered. Panel (c) shows the cooling and indicates strong 
cooling at each of the avoided crossings. 



high phonon occupancy) there is excellent agreement between semiclassical and quantum 
results. At low pressures (near ground state cooling) however, the quantum spectrum 
shows a characteristic asymmetry, such as was observed recently in experiments on 
photonic cavities [12] . 

3.3. Comparison between perturbation theory, semiclassical and quantum results 

The equilibrium variance (and hence the final phonon number) of the mechanical 
oscillator is, 




thus, the final equilibrium temperature of the mechanical oscillator after 
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(b) Quantum (assymmetric sidebands) 




-4-2 2 
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Figure 6. Maps of the displacement noise spectra S xx {uj), showing mode splitting for 
similar parameters to Fig. [3J near the quantum limit, except that here Ai = —1.5 MHz 
is held fixed while A2 (vertical axes) is swept, (a) Shows the semiclassical spectrum 
which is symmetric in frequency uj. (b) Shows the quantum spectrum which is 
asymmetric. In both cases, triple hybridization appears clearly and is seen near 
A2 fa — 1.0MHz. (c) shows the corresponding behaviour for A = k/2, showing that 
the triple peak structure has disappeared. Note that \ogS xx (uj) is plotted. 

optomechanical cooling is k%T eq = l/2mu 2 M (X 2 ). Noting the rescaling (X 2 ) = 
2X ZPF {x 2 ) and setting kpTeq = ((n) + 1/2)Hlum, we can write (x 2 ) = (n) + 1/2. 

Using Eqs. ( IT6|) . ( )20|) and ( !T5|) . we can investigate final equilibrium phonon 
numbers (and the minimum achievable for levitated self-trapped spheres) comparing 
quantum, semiclassical and perturbation theory respectively. In Fig. HI we compare 
the corresponding equilibrium phonon numbers, (ti)qm, (n)sc an d (ti)pt respectively 
for the unusual triple cooling resonance region shown in Fig. |2j Cooling to near the 
ground state (n) ~ is possible for a pressure of order 1CT 6 mBar, even for modest 
driving powers of 2mW and values of A ~ 3 x 10 5 Hz corresponding to spheres of order 
r ~ 100 nm. 

4. Strong coupling regimes: triple mode splitting and bistability 

The multimode, or at least two mode, self-trapping regime may permit new possibilities 
for position sensing and for controlling entanglement between two optical modes 
and the mechanical resonator. Here we investigate regimes where such effects are 
clearly apparent. The implications of the measurement for the accessible range of 
optomechanical coupling strengths suggests that multiple hybridization and bistability 
are quite accessible with reasonable cavity parameters. 

In Fig. [5] we investigate the complex behaviour of the eigenmode frequencies of the 
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(a) Cooling Rates (b) Noise Spectra 




0.0 



Figure 7. Mode mixing and bistability for A = 3k = 6 x 10 5 Hz. We consider a 
relatively low input power of 0.37 mW. (a) plots the optomechanical cooling rate (blue 
indicates cooling, brown indicates heating). The red dashed line indicates the locus 
of bistability as a function of the detunings Ai and A2. The discontinuity in the 
cooling can be discerned near the strong doubly resonant cooling region, (b) shows 
displacement noise spectra S xx (uj) as a function of oj plotted along Ai = — 1 MHz 
(along black horizontal line in upper panel). We sweep in increasing A2. At the point 
where Ai = — 1 MHz intersects the bistability, a discontinuity in the noise spectra is 
apparent. On one side of the discontinuity there is strong hybridization between the 
mechanical mode and optical mode 2; this changes abruptly across the discontinuity 
to hybridization between the mechanical mode and optical mode 1 for larger A2 . This 
may allow for control of entanglement between the modes. The map corresponds to 
near- vacuum conditions so the system is near the quantum ground state in this regime 
(as evidenced by the asymmetric sidebands). 

self-trapped, levitated system. On the left panels (Fig. (a)) we plot the noise spectra 
of the two optical modes, which exibit sidebands near u ~ um since the corresponding 
optical fields are modulated by the motion of the mechanical oscillator. Here we fix one 
detuning (Ai = —1.15 MHz) and look at the behaviour as the other detuning is varied. 
The sidebands are displaced in frequency and split: one effect is simply due to the 
dependence of u>m on Aj (unique to the levitated system); it arises from the calculation 
of the equilibrium fields and frequencies. The other effect is due to normal mode mixing 
(hybridization of light and matter modes) arising from the linearised equations. If 
um ~ Ai ~ A 2 simultaneous hybridization is observed, provided gi^ > k. Figure [5] 
(b) shows that there are several distinct avoided anti- crossings, where the dominant 
character of each eigenmode changes; if two crossings coincide, the spectra show a 
characteristic triple-peak structure (symmetric about u = in the semiclassical regime 
shown here). Panel (c) shows that the corresponding cooling rate is enhanced at each 
avoided level crossing. 



Dynamics of levitated nanospheres: towards the strong coupling regime 16 

In Fig. E] we plot displacement spectra corresponding to Fig. [3J but over a range of 
values of A 2 . A log-scale is used for S xx (u). The triple mixing which can appear when 
two avoided crossings nearly coincide is clearly apparent at A 2 — 1 MHz. 

Static bistability in a cavity of varying length has been seen experimentally 
[3~T] . The potential for generating entanglement has recently been investigated in an 
optomechanical system [32]; however, a relatively high laser power P ~ 50 mW is 
required. 

For the self-trapped systems, the incoherent sum of the optical standing-wave 
potentials cos 2 (kx — 0i) and cos 2 (fcx — 2 ) does not by itself produce a double-well 
structure; nevertheless, as we see below, in combination with optomechanical shifts, 
bistabilities are observed, even for weak driving. Whether a double-well structure 
emerges, or not, is completely independent of the driving power (where P oc E\) and 
can emerge at very low input powers, as we demonstrate below. It is easy to see that 
the levitated particle moves in an effective static potential V(x) where: 



dV(x) 
dx 



hkAEl 



sin 2(kx — <pi) 
](k/2) - iA 1 (x)\' 



+ R 2 



sin 2{kx — 2 ) 
\{K/2)-ik 2 {x)\< 



(22) 



and 



V(x) = HAE 2 



tan 



-i 



Ai(x) 



R tan 



-i 



A 2 to 



Here, note that the shifted detuning Aj 



x 



Aj + Acos 2 (kx 



(23) 

is dependent on 



x, not the equilibrium displacement Xq. This potential admits two stable equilibrium 
points over parameter regimes where A ^> k (in practice, such bistability is observed 
already for A/ k ~ 3). It is evident that the driving power factors out, so does not 
affect the shape of the potential, providing only a scaling factor. We show in Fig. [7] 
that for a high A/ k ratio simultaneous hybridisation and bistability co-exist: we show 
that that for P = 0.37 mW , A = 3k , R = 0.15 and modest photon numbers n\ ~ 10 8 
we can switch discontinuously from hybridisation between the mechanical mode and 
optical mode 1, to hybrization between the mechanical mode and mode 2. We take 
0i = 0, 2 = 7r/4. In the noise spectra, the switch is heralded by a large zero-frequency 
peak in the displacement spectra, which is clearly apparent in Fig. [3 



5. Experiment 

5.1. Current experimental status: Loading protocols and variation of trap frequency 
with radius 

We have built a simple standing wave dipole trap to develop protocols for loading a 
single nanosphere into the trap to confirm that nanospheres with a range of radii around 
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QPD 




nanosphere XYZ A *. 

Figure 8. (Colour online) Schematic diagram of the standing wave trap. It is formed 
from two counter-propagating 1064 nm beams focused inside a vacuum chamber. Light 
at 532 nm enters via one fibre to image the sphere. Images and measurement of the axial 
and transverse position of the trapped nanosphere as a function of time are measured 
by a CCD camera and quadrant cell photodiode (QPD) respectively, through a long 
working distance microscope outside the vacuum chamber. 

100 nm can be trapped. Importantly, we have measured the variation in trap frequency 
with sphere radius so that realistic values of optomechanical coupling strengths, for a 
given nanosphere radius, can be included in our models. In this section we explain how 
the size- dependent modulation function f(r), used in the theory section above to obtain 
A(r) and the optomechanical coupling strengths, was measured. 

A schematic of the standing wave trap used in our experiments is shown in Fig. |HJ 
The standing wave trap consists of two focused beams that counter-propagate and 
overlap near their foci. The two laser beams, derived from the same laser at a wavelength 
of 1064 nm, enter the trapping region via optical fibres. The light exiting the fibres is 
focused using aspheric lenses (Thorlabs C140TME) with a focal length of 1.45 mm and 
numerical aperture 0.55. The power in each trapping beam after it has passed through 
each lens is measured to be 150 ± 10 mW, and the best focused beam waist (radius) 
is theoretically 1.7/im. To optimize the alignment of the trap, we maximize the light 
through-coupled from one fibre into the other. This is accomplished by mounting one 
optical fibre and its aspheric lens on an XYZ flexure stage. The alignment is done inside 
the vacuum chamber at atmospheric pressure. 

A long working distance microscope (Navitar Zoom 6000 system, with up to 45x 
zoom) is used to image the trapped sphere. The image is split into two using a 
beamsplitting plate, with one image directed to a CCD camera for diagnostics and 
the other aligned onto a quadrant cell photodiode (QPD) which measures position 
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fluctuations as a function of time in two orthogonal axes. We define the axial direction 
as that along which the trapping light propagates, and the transverse direction as the 
orthogonal axis in the focal plane of our imaging system. Light at 532 nm is used to 
illuminate the trapped sphere, as the QPD is more sensitive at this wavelength. The 
green beam enters the system via one of the optical fibres, as shown in Fig. EJ and a 
filter is used to stop 1064 nm light reaching the detectors. The power of the 532 nm 
beam is 10 mW. 



Figure 9. An image of a string of 100 nm diameter beads trapped in the standing 
wave trap. A single bead is trapped by continually blocking and unblocking one of the 
trapping beams until only one sphere is trapped. 

Silica (SiC^) nanospheres, manufactured by Microspheres-nanospheres and Bangs 
Laboratories, are introduced into the trapping region at atmospheric pressure via an 
ultrasonic nebulizer (Omron NE-U22). These spheres range in radius from 26 nm to 
510 nm and are suspended in methanol. The nanosphere solution is sonicated using 
an ultrasonic bath for at least an hour before trapping to prevent clumping. Once 
introduced into the trapping region the methanol surrounding the spheres rapidly 
evaporates and the spheres are trapped over many fringes of the standing wave, as 
shown in Fig. [9j As our imaging system does not have single fringe resolution we cannot 
determine if more than one sphere is trapped in a single fringe by this method. However, 
this information can be inferred from the relative intensity of the light scattered from the 
trapped spheres and also by the reduced stability of the particles in the trap when more 
than one particle is trapped. To reduce the number of trapped particles the trapping 
light is briefly blocked and unblocked. This is repeated until a single sphere is visible 
in the trap. At this pressure, where there is a strong damping force from air the sphere 
can be held in the trap indefinitely. 

To measure the trap frequency the air is pumped from the system, and at this 
point no more spheres enter the trap, as without air-damping their velocity is too high. 
The air pressure in the trap is reduced to 5 mbar, so that clear trap frequencies can be 
obtained from the power spectrum of the position fluctuations of the trapped sphere, as 
recorded on the QPD. Example power spectra are shown in Fig. [10l Above 5 mbar the 
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damping of the motion in the trap due to air broadens the peak in the power spectrum 
so that finding an accurate trap frequency is difficult. Below pressures of 5mbar the 
spheres become unstable in the trap and escape. This is most likely due to radiometric 
forces which have been compensated for in other experiments using feedback techniques 
[251 126] . At 5mbar the damping rate due to gas collisions is significantly less than our 
lowest measured trap frequencies, and thus the measured frequency at this pressure is 
a good approximation to the bare trap frequency which would be measured in vacuum 
without damping. 

The angular axial trap frequency for a small polarizable particle in a standing wave 
is u a = 2nf a = W^ga, where the polarizability of a sphere of refractive index n is 

a = 47reor 3 ^^. The maximum intensity in the radial center of each equal intensity 
beam is given by Jo, and k is the magnitude of the wavevector of each beam. The 
sphere has mass m = 4/3irpr 3 , radius r and density p ~ 2000 kgm~ 3 . The transverse 
trap frequency is given by co t = J ' m 8 { "^ 2 , where w is the focused spot size (radius) 
of the two counter-propagating beams. From these expressions the ratio of the trap 
frequencies is given by td a /tu t = kw/V2. 




Figure 10. (Colour online) Power spectra at 5mbar calculated from a measurement 
of the position of a trapped 200.1 nm diameter nanosphere as a function of time, using 
a QPD. (a) The transverse frequency, and (b) an axial frequency. Outlier points are 
due to electronic noise. Red lines show fitted Gaussian functions, from which the trap 
frequencies are extracted. 

The trap frequencies in each axis are determined by fitting measured position 
fluctuation power spectra using ^-uj^yt+u) 2 ^ ' wnere k B is Boltzmann's constant 
and T is the damping rate. The fit to the data is shown in Fig. [TOl with r /27r ~ 2.4 kHz 
for the axial trap frequency. Several sets of data were taken for different spheres of the 
same nominal radius and the measured axial and transverse trap frequencies for each 
size sphere are shown in Fig. [TU The derived trap frequencies for each sphere radius are 
the average over different experiments at each radius, and the errors are the standard 
errors in the mean. The uncertainty in the sphere radius is taken from the information 
supplied by the manufacturer. Two axial (red and green data in Fig. [11]) frequencies 
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and one transverse frequency (blue data points) were measured. When a particle is 
tightly trapped by the optical field only one axial frequency is expected from a single 
sphere in a standing wave. The lower axial frequency (in green in figure [TTj) is always 
observed in the data and this is taken as the true axial frequency. The higher frequency, 
which is often present in the data, may be due to the trapping of two spheres in a single 
anti-node, with the higher frequency occurring due to optical binding, which requires 
further study [27]. The higher axial frequency also changes rapidly with sphere size, 
indicating that it is not the true axial trap frequency, which should be almost constant 
for the small spheres. The presence of a single axial frequency is, we believe indicitave 
of having trapped a single sphere. 

Although we don't know the radial dimensions of the beam within the trap we can 
estimate this value from the ratio of the axial to transverse trap frequencies for small 
spheres. The spot size from this ratio is w = and for ^-=9.8 this gives a spot 

size of w = 2.3 /zm. Since the trap frequency with two overlapping beams of size 2.3 /im 
would be equal to 207 kHz with a power in each beam of 150 mW, and we only measure a 
maximum axial trap frequency of approximately 40 kHz, we conclude that the particles 
are trapped in a standing wave formed where the waist of one beam is much larger. 
If one spot size is 2.3 /zm the other would have to be 15 /im. A plot of the calculated 
axial trapping frequency, found by calculating Maxwell's stress tensor ([28, 29J), is also 
shown in Fig. [HJ Like the experimental data the trapping frequency is constant for 
small spheres and decreases to approximately zero when the particle size is comparable 
to the size of the interference pattern produced by the standing wave. At larger radii the 
force on the particle changes sign and a stable trap is formed in a node of the standing 
wave, as shown for the particle of radius 510 nm. Our measurements confirm that for 
particle radii less that 200 nm the simple dipole model for the nanospheres is adequate 
for modelling the cooling and dynamics of the nanospheres in an optical cavity utilising 
1064 nm radiation. 

Our experiments have shown that optical traps without feedback are currently 
limited to operation at pressures down to a few millibar for all particles that we have 
measured. In addition this limiting pressure did not change by reducing the intensity by 
50%. This radiometric force is due to localized heating of nanosphere and the subsequent 
heating of the surrounding air. At low pressures, when the mean free path is comparable 
to the size of the nanosphere, the radiometric force competes with and eventually 
dominates the dipole force which traps the particle. While feedback techniques have 
been successful [26], decreasing the absorption of the nanospheres is another route to 
minimising radiometric effects. This is feasible since all the spheres we have used in 
this study are not made of optical quality glass but from colloidally grown nanospheres. 
Finally, we have also successfully trapped silica spheres in an ion trap at pressures of 
10 _6 mbar which could be used to load an optical trap formed by a cavity at lower 
pressures where radiometric forces are not significant. 
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Figure 11. (Colour online) Measured trap frequencies as a function of sphere radius. 
Points plotted in green are the axial trap frequency, blue are the transverse trap 
frequency and the red data points are the higher axial trap frequencies which are 
believed to be due to optical binding. The solid black line is a theoretical curve 
derived from a numerical calculation 28 . 

Conclusions 

We have described a study of the dynamics and noise spectra of self-trapped 
levitated optomechanical systems. We have been able to show, by combining 
experimental measurements and theoretical calculations that strong light-matter 
coupling is attainable over a wide range of particle sizes, and that these can be trapped. 
The interdependence of the mechanical and optical mode frequencies, unique to self- 
trapped levitated systems provides a complex and interesting side-band structure, 
including multi-mode mixing and bistabilities which we aim to explore experimentally. 
These conclusions are supported by measurements of trap frequency made in an optical 
standing trap where we have demonstrated a protocol for loading a single nanosphere 
in a single antinode. 
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